Stochastically driven single level quantum dot: a nano-scale finite-time 
thermodynamic machine and its various operational modes 
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We describe a single-level quantum dot in contact with two leads as a nanoscale finite-time 
thermodynamic machine. The dot is driven by an external stochastic force that switches its energy 
between two values. In the isothermal regime, it can operate as a rechargeable battery by generating 
an electric current against the applied bias in response to the stochastic driving, and re-delivering 
work in the reverse cycle. This behavior is reminiscent of the Parrondo paradox. If there is a 
thermal gradient the device can function as a work-generating thermal engine, or as a refrigerator 
that extracts heat from the cold reservoir via the work input of the stochastic driving. The efficiency 
of the machine at maximum power output is investigated for each mode of operation, and universal 
features are identified. 
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Parallel to spectacular developments in bio- and nano- 
technology, there has been great theoretical interest in 
the study of small-scale machines. A well documented 
case is the small-scale Carnot engine, in which the op- 
erational unit is subject to thermal fluctuations [iHj]. 
Of greater biological relevance are machines that convert 
one form of work to another, and yet these have received 
far less attention In this letter we introduce an elec- 
tronic nano-device that allows several modes of opera- 
tion. The device is a single-level quantum dot subject 
to stochastic driving while in contact with two reservoirs 
that may be at different temperatures and chemical po- 
tentials. Its properties can be derived from a stochastic 
thermodynamic description We investigate in ana- 
lytic detail various operational regimes. When operating 
under tight coupling conditions, familiar features are re- 
covered in appropriate limits: Carnot efficiency for re- 
versible operation when the reservoirs are at different 
temperatures, universal features of efficiency at maxi- 
mum power [3|, 0] , and efficiency at maximum power close 
to the Curzon-Ahlborn efficiency Q. When the reser- 
voirs are at the same temperature, the work done on 
the dot by the switching can reverse the "normal" direc- 
tion (from high to low chemical potential) of the current. 
Thus, the engine can be seen as a technologically relevant 
implementation of the Parrondo paradox [9( in that the 
switching can induce an electron flow against the chem- 
ical gradient. When operating under tight coupling, the 
efficiency at maximum power starts from the universal 
value of 1 1 2 close to equilibrium and increases monoton- 
ically to 1 as one moves further into the nonequilibrium 
regime, and can thus be much higher than in the tra- 
ditional implementations of the Parrondo paradox [Toj . 
The same efficiency is observed when the engine works 
in the reverse mode. 



Model and dynamics - We consider a single-level 
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FIG. 1. The model consists of a single level quantum dot. 
A stochastic external force causes the energy of the dot to 
undergo transitions between two values, e u > £d, at random 
times with rates k + and k~ . The dot also exchanges electrons 
with two leads that may be at different chemical potentials 
\il and \xl and temperatures Tl and Tr. 



quantum dot whose energy is stochastically switched be- 
tween an upper and a lower value, £j with j — u, d. The 
upward and downward rates are k + and k~ . The dot is in 
contact with a left and a right lead, v = L, R, at chemical 
potentials \i v and temperatures T v . The transition rates 
of an electron into lead v from the dot, and out of lead v 
to the dot, are given respectively by fc~. — r„j(l — f V j) 
and = T vj fuj. Here f vj = (1+exp {{e s - ^)/T^}) _1 



is the Fermi distribution in lead v, and T 



is the cou- 



pling strength between this lead and the dot in state 
j. The four possible states of the system are denoted 
by {ul, uO, dl, dO}, where jn defines whether the level 
j = u, d is empty or occupied, n = 0, 1. The Markovian 
master equation for the evolution of the state occupation 
probabilities in terms of fe • = k L j + k R j , the total tran- 
sition rate out of the dot (+) or into the dot (-) from 
either lead, is given by 
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The probability currents between the four states are 
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0,1, and Zj^jo - 
ktjPjo ~ k~jPji for j — u, d. It is convenient to in- 
troduce the total current to the dot when it is down 
or up, that is, the sum Zji<-jo = Y^v-^jU-jO' Us- 
ing (P) we easily verify that at steady state there is 
the appropriate balance between the currents, that is, 

I = I u i4-di = ~ AiQ<-do = —lui^uo = Idi<-da- The 
steady state probabilities can be obtained analytically 
(not shown here) and lead to the following expressions 
for the currents: 

k + k~ \ J2v T^ii (^Ldfld + Tfld/fld — ^dfvu) 
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. k + {Tufhd — ^Lufhu ~ ^RufRu) 



•(2) 



k-T d + k+(T u + T d T u ) 
We have introduced the combination T d = T^d 
and similarly for T u . 2^i]_ u0 ^ s obtained from T dl< _ , i/n 
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with the substitutions u <-> d and k + <-> k~ . 

Thermodynamics - The energy current injected in 
the system by the stochastic driving reads 



Text — (Su — £d)%ul<-dl — (^u — £d)T, 



(3) 



while the matter (M) and energy (E) currents entering 
the system from lead v are given by 
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The heat flux from the lead v is 
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It is easy to verify matter and energy conservation in the 
steady state, iff 
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1% . As tive whichever state the dot is in, that is, X^l_^ < 
for j — u. d. However, when the stochastic driving in- 
duces switching between the two states, the resulting net 
current can invert and become positive, with electrons 
flowing against the chemical bias. This remarkable effect 
can be seen as a new version of the Parrondo paradox [j| . 
Thcrmodynamically, the phenomenon corresponds to the 
transformation of one type of work, that which is involved 
in the modulation of the dot energy level, into another, 
the pumping of electrons from low to high chemical po- 
tential. The power at which this process takes place, and 
the corresponding efficiency, are given by: 
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a result, power becomes the sum of two contributions 
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The first is the contribution of the energy flux injected 
by the stochastic driving. The second is the energy flux 
required to bring an electron from the left lead through 
the dot to the right lead. Since at steady state entropy 
production, Si, is minus the entropy flow, that is, Si = 
— YlvQ /T v > 0, we find that entropy production is 
the sum of three force-flux terms, 
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The system reaches equilibrium when entropy produc- 
tion vanish, Si = 0, which implies that all currents in 
the system also vanish. In general, this requires that the 
three thermodynamic forces vanish separately, i.e., that 
£u = £d, t l L = (ar and Tl = Tr. This is however not nec- 
essary when TLu/^Ld — > (disallowing transitions from 
« to L), and TR d /TR U — ¥ (disallowing transitions from 
d to R). This combined limit corresponds to a regime 
of tight coupling, where the transport of matter, energy 
and heat become proportional to each other because, by 
removing the possibility of transitions from u to L and 
from d to R, there is a single pathway connecting the left 
and right leads. As a result entropy production can be 
expressed as Si = XI, i.e. in terms of the "collapsed" 
effective force 
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T^'/sd = a{f L d - Jru), 
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where a = (k+k'T Ld T 'r u ) / '{(k+ + k-){k-T L d + k+T Ru + 
^LdTRu)}- We now see that equilibrium only requires the 
effective force to vanish, X = 0, without the requirement 
that the three thermodynamic forces vanish separately. 

Current rectifier and Parrondo paradox - We re- 
turn to the general expression ([7]) and take the two leads 
to be at the same temperature, T = Tl = Tr but at dif- 
ferent chemical potentials. Wc define A/i = /jr — ^l > 
and Ae = e u — e d > 0. Entropy production then becomes 
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> 0. 



(10) 

In the absence of the external (stochastic) driving, the 
quantum dot remains in its initial state, either in state 
u or in state d. As a result of the direction of the 
chemical gradient, the electronic current will be nega- 
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We carry the analysis further in the tight coupling 
limit, Ijfcf = X, where the efficiency 77 becomes current 
independent, 77 — A/x/Ae. The two forces appearing 
in Eq. 0, Ae/T and A/x/T, collapse into a single one, 
X = (Ae — A/j)/T. Furthermore, from ((5]) we see that 
the condition for pumping electrons against the bias, 
fhd > Iru, is achieved with the Fermi distributions when 
Ae — A/i > 0, that is, when X > 0. The upper bound 
of the efficiency ([TT]) is r\ = 1 (when A/i = Ae) and is 
reached when entropy production vanishes, X — 0. This, 
however implies that the current and hence the power 
vanishes. This reversible operation is not the useful limit 
to consider in practice. Instead, we next consider the 
efficiency at maximum power. 

We start by analyzing the linear response regime near 
equilibrium, where currents are expanded to first or- 
der in the forces, 1 = LnAe — L12A/X and l\^ = 
i2iAe — L22A/X. The Onsager coefficients can be calcu- 
lated analytically, and the reciprocity relation L12 = L21 
can then easily be verified. Maximizing the power out- 
put V with respect to the chemical bias A/x leads to the 
condition A/x* = AeL\ij '(2L22), which in turn yields the 
well-known result for the efficiency at maximum power 
in the linear regime 0, 
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Tight coupling in the linear regime implies the usual re- 
lations among the Onsager coefficients, Ln — L 2 2 = Ln- 
In this case the efficiency at maximum power (|12l) reaches 
its upper bound, 77* = 1/2. 

To go beyond the linear regime, we assume tight cou- 
pling from the outset and introduce the convenient com- 
binations x L = (e d - Hl)/T and x R = (e u - hr)/T, 
in terms of which the Fermi distribution is f(x) = 
(1 + cxpx) -1 . Using © in (TIT]) , power can be rewritten 
as V = a (T(x L - x R ) + Ae) (f(x L ) - f(x R )). Maximiz- 
ing power with respect to xl and xr leads to xr = -ij, 
and xr + sinhx_R = Ae/2T. Using these results and the 
fact that X = xr — xl, the efficiency at maximum power 
is found to be 
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(13) 



Expanding this result close to equilibrium, we find: rf = 
1 /2 + X 2 /96 - X 4 /11520 + O (X 6 ) . Remarkably, the effi- 
ciency at maximum power increases monotonically from 
the linear regime value rf = 1/2 (X 0) to rf = 1, 
reached when X — > 00, cf. Fig. [2j 

The above calculations can be repeated when the en- 
gine operates in reverse, using the difference in chemical 
potential as input work and the modulation of the energy 
level as output. This leads to 
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FIG. 2. Efficiency at maximum power as a function of the 
thermodynamic force X [cf. 1 13] - 



Proceeding along the same lines as before, we find that 
the resulting efficiency at maximum power in the tight 
coupling regime is again given by Eq. (|13[) . Hence both 
modes of operation, forward and backward, can have very 
high efficiency, suggesting a possible technological inter- 
est of the device. 

Refrigerator and heat engine - While the isother- 
mal operation of our engine as described above is its most 
interesting feature, it is instructive to ascertain that it re- 
produces known behavior under more conventional oper- 
ating conditions, namely, when the leads are at different 
temperatures, say Tl > Tr. For this recovery it is suf- 
ficient to consider the simplified case of equal chemical 
potentials, /x = /xj, = (ir. Entropy production now reads 



T R Si = Ael + ricQ {L) > 0, 



(15) 



where 77c = 1 — Tr/Tl is the Carnot efficiency. Our de- 
vice operates as a refrigerator when the external driving 
extracts heat from the cold reservoir, > 0. The 

power output and coefficient of performance 77 of this 
process are given by 
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When functioning as a heat engine, the machine produces 
net work on the stochastic driving process, i.e. l ex t < 0, 
at the cost of a driving heat flow < 0. The power 
output and efficiency of this transformation are given by 
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In the tight coupling limit, the collapsed force (j8]) 
appearing in the entropy production, Si — XI, be- 
comes X = xr — xl with xr = (e u — /x)/Tr and 
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(ed — h)/Tl- Here we have used 1]^' = 1 and 
£dX. The efficiencies (|T6|) and (fTF]) now reduce to 
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FIG. 3. Efficiency at maximum power, rj* , for the thermal 
engine as a function of Carnot efficiency rjc, in comparison 
with the Curzon-Ahlborn efficiency t/ca = 1 — \A — Vc- 



Turning to the regime of maximum power in the tight 
coupling regime, we first discuss the heat engine. In or- 
der to maximize the output power, V = — (e u — Ea)T = 
Tl (xl — xr(1 — t)c))Z(xl,%r), with respect to xl and 
xr, we need to solve OV/Oxl — 0, dTldxR — 0. The 
procedure is identical to that of Ref. Q- The result- 
ing transcendental equation can easily be solved numer- 
ically and leads to the efficiency at maximum power dis- 
played in Fig. [3] We note that rf increases monoton- 
ically when driven out of equilibrium. It is bounded 
from above by r)c, while the Curzon-Ahlborn efficiency 
tjca = 1 — \/l — Vc provides a rather tight lower bound. 
The transcendental equation can also be solved pertur- 
batively for small rjc, 
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We thus recover the universal value r/c/2 in the linear 
regime [7], as well as the factor 1/8 for the quadratic 
coefficient. This latter result thus again supports the 
universality of this value (for systems with a left/right 
symmetry) [3j. 

Finally, we turn to the efficiency at maximum power 
for the refrigerator in the tight coupling regime. The 
output power Q' R ' = — (X + xl)TrI{X,xl) has a lo- 
cal maximum with respect to the thermodynamic force 
X which cannot be found analytically. We can solve 
dQ^ jdX = pcrturbatively as a power series in X by 
expanding x L = ao+a 1 X+a 2 X 2 +a 3 X 3 +a 4 X 4 +0(X 5 ). 
We find ao = a 2 = a 4 = 0, a\ = —2, 03 = 1/3. Inserting 
the result in Eq. (| leads to 
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where = r/^, 1 — 1 is the efficiency in the reversible 
limit. 

Discussion - We have presented a detailed analysis of 
a stochastically driven single-level electronic nano-device. 
When operating as a thermal engine or a refrigerator, 



our model reproduces all the expected results. However, 
of special interest is the isothermal case, where the de- 
vice can be used as a work to work converter. It can be 
seen as a novel implementation of the Parrondo paradox, 
with electrons moving up in chemical potential under the 
influence of the randomly switching energy level. The 
asymmetry in the system is realized via the tight cou- 
pling condition, which implies that each of the energy 
levels of the dot is coupled to a single heat bath. In this 
case, the efficiency at maximum power is very high (up to 
1), suggesting the potential technological importance of 
this mode of operation. Remarkably, the efficiency of the 
device is equally high in the reverse mode, where work 
is extracted from electrons moving down in chemical po- 
tential. Hence, our device can alternate, for example, 
between conversion of chemical energy into electrical en- 
ergy, and vice-versa, thus operating as a highly efficient 
rechargeable battery. It remains to be seen whether the 
simplifications that allow the detailed analysis presented 
here, such as the fully asymmetric coupling, the absence 
of line broadening, and the weak coupling assumption do 
not significantly reduce the device efficiency, and whether 
the technological challenges that the operation of such a 
nano-scale device present, can be overcome. 
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